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Abstract 

We provide a physical analysis of the charging and detection of the first few electrons in a laterally-coupled GaAs/AlGaAs 
quantum dot (LCQD) circuit with integrated quantum point contact (QPC) read-out. Our analysis is based on the numerical 
solution of the Kohn-Sham equation incorporated into a three-dimensional self-consistent scheme for simulating the quantum 
device. Electronic states and eigenenergy spectra reflecting the particular LCQD conflnement shape are obtained as a function 
of external gate voltages. We also derive the stability diagram for the first few electrons in the device, and obtain excellent 
agreement with experimental data. 

PACS numbers: 73.21.-b, 72.20.My, 73.40.Gk 
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I. INTRODUCTION 

Lateral GaAs/AlGaAs quantum dots (QD's) are now 
routinely fabricated with planar technology. ^Ij Three- 
dimensional (3D) quantum confinement is achieved, in 
part, by using the GaAs/AlGaAs semiconductor het- 
erostructures to confine the conduction electrons into a 
two dimensional electron gas (2DEG) at the interface be- 
tween the two materials. By placing metal gates on top 
of such a structure, carrier confinement in other in-plane 
directions can be realized by energizing the gates that 
create lateral energy barriers to electrons in the 2DEG. 
Design of these QD's, which previously contained tens 
of electrons, has been improved to operate them in a 
few-electron regime where the charging of the very first 
electrons can be observed experimentally. Two quan- 
tum dots can be placed adjacent to each other to form 
a laterally coupled device with both electrostatic and 
quantum- mechanical coupling between them. ^ Fine 
variations of the top gate biases change the confinement 
of each dot, while precise coupling between them through 
the central gates leads to a fully tunable two-qubit quan- 
tum system, which can be used as a building block for 
quantum computing. 5] 

Recently, it has been shown that laterally-coupled 
quantum dots (LCQD) containing a few conduction elec- 
trons could be coupled to single charge detectors to form 
an integrated quantum circuit. Q The read-out of the 
charge state in the LCQD is realized by integrating mono- 
lithically quantum point contacts (QPC's) adjacent to 
each of the QD's. Each QPC can be calibrated through 
electrostatic coupling with the dots so that its conductiv- 



ity changes abruptly once a single-electron charging event 
occurs in one of the dots.|^|3 With this sensitive detec- 
tor, it is then possible to obtain the "stability diagram" 
that describes the stable charge regimes of the LCQD as 
a function of the tuning (plunger) gate biases. ^3, '5] 

This quantum dot circuit has a two-fold advantage: it 
is possible (i) to scale it to a quantum dot array, (ii) to 
perform single-quantum sensitivity measurements, both 
of which are favorable features of a realizable quantum 
computer. 

In this paper we study the properties of the above 
circuit via numerical simulation that involves the self- 
consistent solution of coupled Poisson and Kohn-Sham 
equations discretized on a 3D mesh.[TlllT^ In Sec. II we 
describe the LCQD structure and in Sec. Ill we present 
the approach for solving the Kohn-Sham equations in the 
device environment within the local spin density approx- 
imation (LSDA) and express the criterion used to de- 
termine the charging events as a function of the applied 
gate biases. In Sec. IV we present our simulation results 
of the circuit, including both electrostatic and quantum- 
mechanical features, the functionality of the QPC's, and 
the stability diagram in the few-electron charging regime. 
Finally, we summarize our work in Sec. V. 



II. DOT STRUCTURES 

Figure 1(a) shows the top view of the LCQD and QPC 
gates in the a;y-plane. Q Top L-, R-, T- and M-gates 
are used to define the two coupled-dot region. Among 
them, the T- and M-gates can also control the cou- 



pling between the two dots. The PL- and PR-gates, 
called the "plungers," have smaller feature sizes than 
the other gates and are used for fine tuning the con- 
finement of each dot. The QPC-L and QPC-R gates 
are associated with the L- and R-gates (via the tips) 
to form the QPC detectors. Charging paths into the 
dots (shown by the ovals) from external reservoirs are 
shown by curved arrows, whereas the QPC currents are 
shown by straight arrows. Figure 1(b) shows a cross- 
sectional view of the layer structure in the z-direction. 
Our model involves four different layers of semiconduc- 
tor materials (from top to bottom): a 50-A thick n-typc 
{Nd = 1.5 X 10^* cm~3) QaAs layer, a 650-A thick n- 
type {Nd = 0.31 x 10^* cm"^) Alo.27Gao.73As layer, a 
200-A thick undoped Alo.27Gao.73As layer, and a 1610- 
nm thick p-type (Na — 1.0 x 10^^ cm~'^) GaAs layer. 
The 2DEG is formed at the interface between the un- 
doped AlGaAs layer and the lightly p-type doped GaAs 
layer (900 A below the top surface). 



III. NUMERICAL MODEL 

The electron density in the LCQD region is obtained 
by describing the charge carriers within the density func- 



tional theory (DFT) that incorporates many-body effects 
among particles. [13I In order to take into account the 
spin dependence of the electron-electron interaction, the 
Kohn-Sham equations 14] for spin-up (f) and spin-down 
(I) are solved simultaneously: 



i/Vj(r)-4v4(r), 
i?Vi(r)=£fv4(r). 



(1) 



Here and "0^ are the corresponding eigenenergies 
and eigenfunctions of the Hamiltonian H^^^'>: 



m*(r) 



ir) + AE, + (j,l^P{n), (2) 



where m*(r) is the position dependent effective mass. 
(/)(r) — 4>ext+4'ion-\-4>H is the clectrostatic potential which 
consists three parts: (j^ext is the potential due to external 
gate biases, is the potential resulting from ionized 
donors and acceptors, and (j)H is the Hartree potential 
accounting for repulsive electron-electron interactions. 
lS.Ec is the conduction-band offset between different ma- 
terials, and <j3'x}\v) is the exchange-correlation poten- 
tial energy for spin-up (t) and spin-down (|) computed 
within the local spin density approximation (LSDA) ac- 
cording to Perdew and Wang's formulation. |15j Hence 



FIG. 1: (a) Layout of the top gates (Light gray areas show 
the gate pattern for the LCQD and the QPC's; ovals show the 
dots; curved arrows show the possible charging current paths; 
and straight arrows show the QPC currents.), (b) Layers of 
the heterostructure (not to scale), after Elzerman et al.. 4] 



our approach is spin unrestricted by allowing for differ- 
ent orbitals with different spins. 

The electron density rt(r) in the LCQD region is 

|2 I 2 

1=1 i=i 

where iV| iV^ = iV is the total number of electrons in 
the dots. 

The electrostatic potential (f){v) is computed by solving 
Poisson's equation 

V [e(r)V0(r)] = -p(r), (4) 

where e(r) is the position-dependent permittivity and 
p(r) is the total charge density given by 

p{v)=q [N+ (r) - (r) + p(r) - n(r)] . (5) 

Here A'^^(r) and N^{r) are the ionized donor and ac- 
ceptor concentrations in the relevant device layers, p(r) 
is the hole concentration, and rt(r) is the total electron 
concentration given by Eq. (3) in the QD region, while 
outside this region the free electron charge is entirely 
determined by using the semi-classical Thomas-Fermi 
approximation, 

We solve Kohn-Sham and Poisson equations self- 
consistently by finite element method. 11, 12] Zero nor- 
mal electric field on lateral and bottom surfaces and 
Schottky barrier values on the top surface are imposed 
as boundary conditions for the solution of Poisson equa- 
tion. Since the quantum dots are much smaller than the 
physical dimensions of the device, the wavefunctions ac- 
tually vanish long before reaching the device boundaries. 
This allows us to embed a local region in the global mesh 
for solving the Kohn-Sham equations. This local region 
is chosen large enough to ensure vanishing wavefunctions 
on its boundaries. A non-uniform 3D grid of 141, 52 and 
71 mesh points in the x-, y- and z-directions respectively, 
is used for solving Poisson equation, while 71 x 45 x 19 
grid points are used to discretize the local region where 
Kohn-Sham wavefunctions are evaluated. 

Because the LCQD are weakly coupled to the external 
reservoirs, we assume that electrons in the dots are com- 
pletely localized in that region. At equilibrium, and for 
a given bias, an integer number of electrons N minimizes 
the total energy Et of the dots. In order to determine 
N , we use the Slater formula: [l^ 

Et{N+1)-Et{N)^ [ eLuo{n)dn ^ eLuo{l/'2)-EF, 

JQ 

(6) 

where Et{N+1), Et{N) are the total energies for iV -1-1, 
N electrons in the dots, and ELUoi^/'^) is the eigenenergy 
of "the lowest unoccupied orbital" with half occupancy. 
The sign change of the right-hand side of Eq. (6), as a 
function of the tuning gate voltage, determines the elec- 
tron occupation in the LCQD. In our simulation, we use 
a variation of the above rule where charging occurs when 
eLf7o(l) — Ep = Ep — eL(7o(0), which was justified in 
Ref. 12. 



IV. RESULTS AND DISCUSSIONS 

Figures 2 shows the conduction band edge profiles in 
the xy-plane at the 2DEG interface (contour plot, Fig. 
2(a)) and in the z-direction (Fig. 2(b)) under the con- 
dition Vl = Vr = Vqpc-l = Vqpc-r = Vm = -0.585 

V, Vt ^ -0.9 V, VpL = VpR ^ -0.15 V (These voltages 
correspond to point A in Fig. 8.) and zero electrons in 
the dots. The Fermi level is set at zero throughout the 
device at the temperature T = 4 K. The LCQD region 
and the QPC region with low equipotential-line density 
are clearly visible in Fig. 2(a). The outer energy barrier 
for the LCQD is ~ 110 meV whereas the energy barrier 
between the dots is ~ 9 meV. A large negative T gate bias 
is used to prevent the wavefunctions from leaking into the 
external reservoirs, which clearly defines the LCQD re- 
gion. The confinement along the z-direction is achieved 
by a quasi-triangular shaped well shown in Fig. 2(b), for 
which the relaxation of the potential to zero-field is not 
shown at the far end (substrate) of the device. Due to 
the strong confinement in the triangular well, only the 
ground state along the z-direction is occupied (the shape 
of the ground state wavefunction along the z-direction is 
shown in the inset in Fig. 2(b).). Under the above con- 
dition, the wavefunction contour plots in the xy-plane at 
the 2DEG interface are shown in ascending energies for 
the first eight spin-up (t) eigenstates in Fig. 3. A sim- 
ilar set of wavefunctions is obtained for the spin-down 
(1) eigenstates (not shown). They are similar to orbitals 
observed in diatomic molecules: the two columns repre- 
sent the familiar bonding and anti-bonding state pairs. 
Notice that the shape of the wavefunctions reflects the 
shape of the confinement seen in the local minima of the 
conduction band edge in Fig. 2(a). 

In Figure 4(a), we show the variation of the first eight 
spin-up (I) eigenenergies when the plunger gate bias con- 
figuration is changed from the values Vp^ = VpR — 
-0.15 V to the new values Vpl = -0.15 V, Vpr = -0.06 
V. The first eight eigenenergies are separated into two 
groups, one for the right dot (solid lines) and one for 
the left dot (dashed lines), which arc lowered simultane- 
ously as the right plunger gate bias increases. However, 
the eigenenergies of the right dot decrease more rapidly 
than those of the left dot because of the proximity of 
the former to the varying plunger. At Vrr — —0.074 
V, the charging of the first electron (spin-up (t)) occurs 
in the right dot, which is indicated by a discontinuity of 
8.2 X 10"'' eV in the variation of the ground state energy 
level with respect to the right plunger gate bias. At the 
same gate bias, we also observe a jump of the conduc- 



FIG. 2: Conduction band edge profile in the LCQD-QPC 
structure (a) contour plot in the xy-plane at the 2DEG inter- 
face (The dashed rectangle shows the location of the dots.), 
(b) along the z-direction with the inset showing the shape of 
the ground state wavefunction {Vpl = Vpr — —0.15 V, zero 
electrons in the LCQD). 



FIG. 3: Contour plot of the first eight spin-up (|) eigenstates 
in ascending energies in the xy-p\ane at the 2DEG interface 
with zero electrons in the LCQD {Vpl = Vpr — —0.15 V). 
The sy-coordinates are given for the lower left wavefunction, 
which is a zoom-in region corresponding to the dashed rectan- 
gular region in Fig. 2(a); all the other wavefunction contour 
plots in this paper are on the same scale. 



FIG. 4: (a) Eigenenergy spectrum (spin-up (|) states) as 
a function of the right plunger gate bias (solid lines: right 
dot; dashed lines: left dot), a, f3 and 7 are three "anti- 
crossing" points, (b) Variation of the conduction band edge 
in the constriction of the left and right QPC's as a function of 
the right plunger gate bias from point A to B in Fig. 8 {Vpl 
is fixed to —0.15 V; the vertical axis of Fig. 4(b) is shifted up 
by 0.0201 eV.). 



tion band edge in the constriction of the two QPC's, i.e., 
2.6 X 10"^ eV for the left QPC and 5.4 x 10"^ eV for 
the right QPC (see Fig. 4(b), where the vertical axis is 
shifted up by 0.0201 eV for clarity). The up-shift of the 
conduction band edge in the QPC constriction results 
from the Coulomb interaction between the electrons in 
the LCQD and electrons in the QPC's, which reduces 
the total charge number in the conduction channel and 
leads to a discontinuity in the QPC current observed in 
experiments. Uj Obviously, the right QPC is more sensi- 
tive to the single-electron charging because of its prox- 
imity with the right dot. On the stability diagram (Fig. 
8), this transition is represented by the vertical A to B 
line with the diamond indicating the charging point for 
the first electron. 

From the eigenenergies variation vs. Vpr diagram 
(Fig. 4(a)), we also observe three "anti-crossing" points 
between the two different sets of eigenenergy levels, each 
arising from the distinct QD's as mentioned above and 
indicated by arrows in Fig. 4(a), i.e., (a) at Vpp ~ 
—0.09 V, between the 3rd and 4th excited states; (/?) 
at VpR ~ —0.07 V, between the 3rd and 4th excited 
states; and (7) at Vpp ^ —0.11 V, between the 5th and 
6th excited states. The behavior of the system near the 
"anti-crossing" points can be further illustrated by ex- 
amining the evolution of the wavefunctions for the "anti- 
crossing" levels. "Interchange" of the wavefunctions is 
clearly observed before and after these points. In Fig. 5, 
contour plots of the wavefunctions in the xy-plane at the 
2DEG interface are shown for the three "anti-crossing" 
points: the 3rd and 4th excited states labeled ai, a2 at 
Vpj^ = -0.10 V and a[, a'^ at Vrr = -0.08 V, respec- 
tively; the 3rd and 4th excited states labeled (32 at 
Vpj^ = -0.074 V and I3[, at Vrr = -0.06 V; and the 
5th and 6th excited states labeled 71, 72 at Vrr, = —0.12 
V and 7^, 7^ at Vrr = -0.10 V. 

The detection of single-electron charging events can 
also be carried out for the B to C transition in Fig. 8, 
in which the right plunger gate bias Vrr is fixed to be 
—0.06 V while the left plunger gate bias is changed from 



VpL = -0.15 V to VpL = -0.06 V. The variation of 
the spin-up (|) eigenenergies with respect to the Fermi 
level and the conduction band edge in the constriction 
of the two QPC's are shown in Figs. 6(a) and (b), re- 
spectively. In this case, the transition of the charging 
state is from one electron in the right dot to two elec- 
trons, one in each dot occupying an individual IS'-likc 
orbital, 10 and occurs when the left plunger gate bias 
VpL is at —0.097 V. The charging of the second electron 
into the system is localized in the left dot and is indicated 
by the jump of the first excited state energy level. Note 
that in this case, the variation of eigenenergies in the left 
QD (dashed lines) is larger than those in the right dot 
(solid lines) . In our LSDA approach, the second electron 
has the same spin (spin-up (|)) as the first one as they 
are uncorrelated by the height of the coupling barrier. 
The corresponding jump of the conduction band edge is 
5.6 X 10-6 eV for the left QPC and 2.8 x lO'^ eV for the 
right one. The left QPC is more sensitive to the second 
electron charging because it occurs in the left dot. Fig. 
7 shows the spin- up (t) wavefunctions in ascending ener- 
gies (from top to bottom) after the charging of the second 
electron for two bias conditions, 1) VpL = —0.08 V and 
2) VpL = -0.06 V (for both cases, Vpr = -0.06 V). 
It demonstrates the evolution of the wavefunctions from 
an asymmetric configuration with different eigenenergy 
levels for two electrons in the system to a symmetric one 
where eigenenergy levels are fully degenerate. 

Following the same procedure as described above, we 
can find another charging path for the first electron 
charging, i.e., from point E to F in Fig. 8, and for the 
charging from one to two electrons, F to G in Fig. 8, 
for distinct stable charge regimes of electrons in the two 
dots. On the path E to F {Vpl is fixed to -0.125 V, Vpp, 
is changed from —0.125 V to —0.07 V), charging hap- 
pens for the first electron (spin-up (t)) in the right dot 
at VpR = -0.082 V; on the path F to G {Vpr is fixed 
to -0.07 V, VpL is changed from -0.125 V to -0.07 V), 
charging happens for the second electron (spin-up {])) in 
the left dot at Fpl = -0.092 V. 

We can further interchange the plunger gate biases and 
obtain different transitions, i.e., from A to D to C and E 
to H to G, as shown in Fig. 8 to realize closed cycles of 
charging and discharging paths. These two closed paths 
(dashed and dotted lines) are shown in Fig. 8. Each 
corner of the two squares is in a different stable charge 
state with numbers in the parentheses showing the elec- 
tron number in the left and right dots, respectively, e.g., 
(0,1) means zero electron in the left dot and one in the 
right dot. On each path, we record the charging points 
(diamonds in Fig. 8) and make linear extrapolations be- 
tween the two charging points on each two parallel paths. 



FIG. 5: Wavefunction (for spin-up (|) states) "interchanges" 
at the "anti-crossing" points corresponding to (a) point a, (b) 
point P and (c) point 7 in Fig. 4(a). 



FIG. 6: (a) Eigenenergy spectrum (spin-up (|) states) (solid 
lines: right dot; dashed lines: left dot) and (b) variation of 
the conduction band edge in the constriction of the left and 
right QPC's as a function of the left plunger gate bias from 
point B to C in Fig. 8 {Vpr is fixed to —0.06 V; the vertical 
axis of Fig. 6(b) is shifted up by 0.0201 eV.). 



FIG. 7: Evolution of the first eight spin-up (j) wavefunctions 
from bias condition 1) Vpl = -0.08 V to 2) Vpl = -0.06 V 
(for both cases, Vpr = -0.06 V). 



which leads to four lines crossing at two points (circles in 
Fig. 8). 

The two crossing points are linked afterwards. Now, 
five segments (solid lines in Fig. 8) separate the dia- 
gram into four regions to define the stability diagram for 
the LCQD system in the few-electron charging regime. 
Each region, separated by the solid lines, indicates a sta- 
ble charge configuration assumed by the LCQD under a 
particular range of plunger gate biases. More interesting 
are the two crossing points (circles), called the double- 
triple point, 0] occurring at Vpl = Vrr = —0.0924 V 
for the three charging states (0,0), (0,1) and (1,0) and 
at VpL = VpR = -0.0847 V for (0,1), (1,0) and (1,1) 
states. We then determine the voltage range of the right 
plunger that spans the distance between the double-triple 
point to be AVp^j = 7.7 mV, which is comparable to the 
experimental result ^ 7.4 mV.Q 

Finally, from the charging diagrams in the few-electron 
regime we extract the addition energy for the second elec- 
tron charging in the right dot: we determine the Vpp- 
voltage interval on the stability diagram for the (0,1) 
configuration {i.e., between the (0,0) configuration and 
the (0,2) configuration in the singlet state) to be 0.1 V, 
which is in excellent agreement with the experimental re- 
sult ^ 0.1 V. 4] By linear projection of this Vpp interval 
to the energy scale, 12j we obtain then the addition en- 
ergy for charging the second electron, which is 2.5 meV. 
By comparing this value to the experimental result of 3.7 
meV,|3 we attribute the difference to the fact that our 
simulation is performed on a coupled-dot system, while 
the experimental result is obtained by grounding one of 
the dots where the confinement is stronger in an individ- 
ual dot compared to our simulation case. 



V. CONCLUSION 

We performed numerical simulations of the electro- 
static and quantum-mechanical characteristics of a novel 
laterally-coupled quantum dot circuit with integrated 
quantum point contact read-out. We were able to repro- 
duce detailed single-electron charging behavior of the el- 
ementary quantum circuit and demonstrate the function- 
ality of the QPC's as single-electron charging detectors. 
In particular, we obtained excellent agreement with the 



FIG. 8: Stability diagram for the first two charging electrons 
characterizing the double-triple point (shown by circles). 
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